Statistical dynamics of a non-Abelian anyonic quantum walk 
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We study the single particle dynamics of a mobile non-Abelian anyon hopping around many pinned anyons 
on a surface. The dynamics is modelled by a discrete time quantum walk and the spatial degree of freedom 
of the mobile anyon becomes entangled with the fusion degrees of freedom of the collective system. Each 
quantum trajectory makes a closed braid on the world lines of the particles establishing a direct connection 
between statistical dynamics and quantum link invariants. We find that asymptotically a mobile Ising anyon 
becomes so entangled with its environment that its statistical dynamics reduces to a classical random walk with 
linear dispersion in contrast to particles with Abelian statistics which have quadratic dispersion. 
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Anyons are point like particles with more general statis- 
tics than bosons or fermions. They were shown to exist in 
systems where the physics is constrained to two dimensions 
[1]. Beyond mere possible existence they where found to be 
a good description for low lying quasi-particle excitations of 
fractional quantum Hall systems [2, 3] and they exactly de- 
scribe excitations in various strongly correlated two dimen- 
sional spin lattice models [ , ]. Recently there has been 
tremendous experimental progress in preparation and control 
of systems capable of exhibiting topological order [ , ] with 
the goal to observe anyonic statistics. This is further motivated 
by the discovery that braiding some types of non-Abelian 
anyons can be used for naturally fault tolerant quantum com- 
puting [ ]. The quantum physics of anyonic systems is very 
rich but is only beginning to be explored in its own right. 
For example, there have been investigations of the equilibrium 
properties of dynamically interacting, but static, non-Abelian 
anyons in chains [ ] and two dimensional lattices [10, 11]. 

Here we describe a simple model which captures some of 
the non-equilibrium physics of moving non-Abelian anyons 
interacting purely due to particle statistics. In Ref. [12] the 
authors introduced a general protocol for quantum walks with 
anyons to describe the dynamics of one mobile anyon braid- 
ing about other pinned anyons on a surface. It was shown that 
while for Abelian anyons the dispersion is quadratic as in the 
usual quantum walk, the non-Abelian walk appears to have 
richer behavior. Transition from coherent quantum to classical 
random behavior with linear dispersion occurs when a quan- 
tum walk strongly decoheres due to interaction with an envi- 
ronment [ ]. Indeed, Ref. [14] found that if one introduces a 
new coin at every other step, or fewer, in a quantum walk the 
dispersion is quadratic but if a new coin is introduced every 
time step then the dispersion is linear. One might expect that 
statistical interactions of non-Abelian anyons would be suffi- 
cient to induce such a transition. We show that this is the case 
by expressing the statistical dynamics of the mobile anyon as 
a function of topological invariants of the links of the anyonic 
worldlines generated during the quantum walk. Specifically, 
we solve for the asymptotic distribution of the Ising model 



non-Abelian anyons, a, which appear as quasi-particle excita- 
tions in the Pfaffian wave function description of the 5/2 filled 
fractional quantum Hall state [3, 15]. It is the purpose of this 
work to determine the behaviour of such a system by analyt- 
ical methods, thus opening the way for modelling complex 
systems that are of interest to statistical physics [16]. 

The setup (see Fig. la) is a surface with n vacuum pairs 
of anyons of topological charge a, and one member of each 
pair participates in the dynamics. The other n anyons are 
moved out of the way or could be excitations on the boundary 
[26]. The participating anyons are canonically ordered on the 
surface with the n—\ pinned anyons and one mobile walker 
anyon. The walker hops between neighboring sites with spa- 
tial index s=l,...n—l and has an additional spin- 1/2 degree 
of freedom (DOF) we dub a coin. The pinned anyons are lo- 
cated in between the sites. The total Hilbert space decomposes 
as H = i^pace <8> H:oin <8> fusion ^ C n ~ 1 <g> C 2 <g> C D which be- 
comes infinite dimensional in the asymptotic limit. Fusion 
DOFs denumerate the number of distinct measurement out- 
comes of topological charge when pairs of anyons are fused 
together [ ]. Its size grows like D~d n where d is the quantum 
dimension of the anyons. For n Ising anyons with total trivial 
charge, we have d = y/2 and D = 2 n / 2_1 , thus effectively they 
introduce half a new coin DOF per step of the walk. 

The dynamics is modelled by a composition of two discrete 
unitary steps W = TU where U acts on the coin and T is a 
conditional braiding operator. It moves the walker to the right 
or left depending on the coin state: 

T = I^l 2 k-l) S pace^l^|0)coin(0|^Vl 
+ k+ 1 )space(^l^l 1 )coin( 1 l^^ 

where {b s }"~l is a set of unitary generators of the braid group 
*B n . The particular representation depends on the braiding and 
recoupling rules for the chosen anyons. Notice that the chiral- 
ity of the mobile anyonic charge current is fixed counterclock- 
wise by this walk. To make T unitary we assume periodic 
boundary conditions (|0) space = \n— l) space ) but will be con- 
cerned with walks satisfying \n/2\ < t so that winding around 
the surface is not an issue. 
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The system's 



K = TlDspace l°)coin I 



initial state is |¥(0)) = 
Infusion wnere \&) i s me vacuum 
configuration of the n pairs of anyons with half the mem- 
bers braided to the right. After t iterations, the state is 
= W f |¥(0)) and the reduced state of the spatial DOF 
of the walker is 



PspaceW = tr coin trf usion 

= /£tr^tr^|2|a|-r + jo)(2|a / |-r + j |(l) 



where 



= (m=i Pa r U) |0) coin (0| (FlU Pjtf) and JL, = 
Bt\&) (<&\Bl,^ . where |<3>) corresponds to the Markov trace 
state (see Fig. lb). The coin histories are given by the vectors 
a,a! G {0, 1}®*, and the projectors for each outcome are P a . = 
I* 2 ./) coin ( a J I ' ^ ne ^ ra ^ wor d f° r a given coin history is 



t-\ 



so+at- 



r+2a?pr 1 «y) 



For a Hadamard coin flip operation £/ 



V 



1 -ly 



(2) 



and 



tr ^ = ^(-l)^ where z(2,a') = Zfr 
is the sum of pairs of consecutive right moves (or 1 outcomes 
of the coin). The trace over the fusion DOF can be related 
to the Kauffman bracket of a link, denoted (L) for a link L. 
Here we get a link which is the Markov trace over the braid 
words for the forward and backward time evolution histories 
as dictated by Eq. 1 : 



Markov 



(3) 



Moreover, 



Markov 



where d is the quantum dimension of the anyons [12]. 

Concerning ourselves with the diagonal ele- 
ments of the spatial probability distribution p(s,t) = 
(2s — t + so | P space (0 \2s — t + so), we have the constraint that 
\a\ = \a!\, i.e. the final position of the walker for the braids 
B$f is the same, and the trace over the coin DOF is non 
zero only if a t = a! t , i.e. the final step is in the same direction 
for both paths. The result is: 



p(s,t) = 



1 



{a,a f ; \a\ 



E 

=1*1= 



(L) (4) 



which for d = 1 and (L) = 1 reduces to the usual quantum 
walk distribution. 

The problem of computing the distribution thus reduces to 
computing statistics of a quantum link invariant. The Kauff- 
man bracket of a link (L) (A) is a Laurent polynomial in the 
argument A that is an invariant for framed, unoriented links 
and is framing dependent. It can be related to the Jones poly- 
nomial Vl, which is an invariant for framed, oriented links but 
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FIG. 1: Anyonic quantum walk in a system with n vacuum pairs 
of anyons in state |<I>). (a) Half of each vacuum pair (depicted as 
connected by a string) participates with one walker anyon carrying 
a spin braiding around n—\ pinned anyons. (b) Link representation 
of the world lines for a Markov closed quantum trajectory that con- 
tributes to the spatial distribution p(3, 5). Here a = (1 , 0, 0, 1 , l) T and 
a' = (0, 1 , 1 , 0, 1 ) T . The link shown is proper and has one Borromean 
ring, i.e. if any of three linked components were cut the others would 
become disentangled. Here arf (L) = — 1 . 



(a) 



FIG. 2: (a) The writhe of a link is the difference of the number of 
positive crossing and negative crossings: w(L) = #£+ —#£-. (b) The 
linking number of two components, shown here as thick and thin 
components, is determined by the relative number of four types of 
crossings lk{L thin ,L thick ) = Mi+Mii^Ma = #£ x _#£ 4 = #£ 2 -#£ 3 . 



which is framing independent. The relation is achieved by in- 
troducing an orientation to each component in L and suitably 
normalizing: 

(L)(A)\ A ^. 1/4 = (- q y 4 r iL) VL(q) 

in such a way that the framing dependence of the Kauffman 
bracket is exactly canceled by the multiplicative factor involv- 
ing the writhe w(L) of L (defined in Fig. 2a). The Jones poly- 
nomial of a link L with variable q = e l2n ^ ^ +2 ) was shown by 
Witten [ ] to be equal to the expectation value of the product 
of the path ordered Wilson loops along the components of the 
links of L in SU (2)& Chern-Simons theory. 

At the special value q = i(k = 2), where the links represent 
braiding of anyons in the Ising model, the Jones polynomial 
of a link can be related to a simpler knot invariant known as 
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the arf invariant [18]. Specifically, 

Vl(0 = 



^(#(L)-l) ( _ 1)arf(i) ifLproper 

if L not proper 



where the number of components #(L) = n here. An oriented 
link is proper if each component L k evenly links the union 
of other components, i.e. ^^/^(L^Lfc) = mod 2Vj. The 
linking number of two components is defined in Fig. 2b and 
can be computed in polynomial time. The advantage of this 
expression is that arf (L) G {0, 1} of a link can be computed 
in polynomial time in the crossing number of a braid presen- 
tation. Hence, unlike the generic case, the Jones polynomials 
at value q = i can be evaluated in polynomial time [ ] . How- 
ever, the number of links contributing to the weight p{s,t) is 

C^ 1 ) + C-i) wri ich is exponential in t so an efficient com- 
putation is not a priori available. 

In the special case the link L has all pairwise linking num- 
bers even, i.e. it is totally proper, then there is a three local 
formula for the arf invariant: 

arf(L)=£ci(L f ) + £c 2 (Lf,L 7 -)+ £ c 3 (Li,Lj,L k ) mod 2 

i i< j i< j<k 

(5) 

where c s (T) is the coefficient of z s+l in the Alexander- 
Conway polynomial of the s component sublink T [18, 20]. 

There is structure to the links in the anyonic quantum walk 
trajectories that simplifies Eq. 5. First, if L is proper than it 
is totally proper. This is due to the following: if we consider 
the set of all the sums S = {£ k ^jlk(Lj,L k )} n j =v this is equal 
to S = {lk(Lj,L w )}j-L w U {£ k ^ w lk(L Wl L k )} where L w is the 
walker's component. The condition of L being proper is that 
every member of S is an even integer. Now, the total link- 
ing number is Y<i<jlk( L ii L j)) = Hk^w l K L w,L k )) = 0. This 
quantity is zero since for the forward half of the braid, i.e. Bl, 
all contributions to links are positive (counterclockwise braid- 
ing) and for the latter half, Bt,\ there is the same number 
of clockwise braids all of whom contribute with a negative 
sign, so S = {lk(Lj : L w )}j^ w U {0}. But the condition that 
L is totally proper is precisely that every member of the set 
{Ik^Lj^Lw)}^ is an even integer. Hence if L is proper then 
it is totally proper. Second, the writhe of any link L is zero 
since the total linking number is zero and the braids act on 
components all with the same orientation (see Fig. lb). 

With regard to the terms in Eq. (5), note that due to causal- 
ity of the worldlines there is no self linking, i.e. c\ = 0. The 
pairwise contribution to the arf is also zero since the sum of 
C2(Lj,Lk) is even as we now show. It is only necessary to con- 
sider links involving the walker L w and every other component 
Lj since the non-walker components are not directly linked. 
Such a link is the braid closure of a braid in the two com- 
ponent braid group (B2 with one generator b so that link can 
be written, (L W ,L/) = (b m ) Markow where m = 2 x lk(L w ,Lj). 
To compute the two point invariant we use the defining Skein 
relation for the Alexander-Conway polynomial 



with the notation defined in Fig. 2a. Moreover, = 1 which 
states that if one of the components is unlinked from the others 
then it can be removed with a multiplicative factor of 1 . The 
polynomial for any link L can be written as V^(z) = YT =0 aiZ l . 
For our pairwise component links (L Wj Lj) one can solve for 
the polynomial by recursion: 



V (L VV ,L J -) =f\m\(z)+f\m\-l(z)> 



where 



fm(z) 



2 m+l ^/ z 2 +4 



((z-v / ?T4)((z + V / ?T4r 



From the coefficient of the cubic term of V( L ^ L ^ we obtain 

c 2 (L w ,Lj) = lk(L w ,Lj)(lk(L w ,Lj) 2 - l)/6. 

Note that the linking numbers are always even, lk(L Wl Lj) = 
2rij with rij an integer, so the sum can be written as 

^ j<k c 2 (Lj,L k ) = §I# w ^-§E y y w n/ 
= 8E^w(% +1 ) ^2N. 

Here we used the facts that the sum of the linking numbers 
is zero, so that Y.j^w n j ~ 0' an d the binomial coefficient is 
always an integer. Finally, the triple component invariant 
C3(L rj L Sj L t ) is known as the Milnor invariant. It counts the 
number of Borromean rings in three component sublinks [21] 
and has nontrivial contribution to the arf of our links. 

This analysis allows us to express the probability distribu- 
tion in terms of simple properties of the anyonic walk: 



p(s,t) = 



{H/e{0,l}^; 
\a\ = \a f \=s,a t =a f t } 



(_ l \z(a,a f )+x(a,a f ) 

- — L — jt ^ proper 

L not proper 



(7) 



(6) 



where the link L is defined in Eq. (3) and the sum of Mil- 
nor invariants over sublinks (L n L s ,L t ) of L is x(a,a 7 ) = 
Hr<s<t c 3(Lr,L s ,L t ). When all links are proper and T(a,a / ) 
is even then p(s,t) —> Pqw(s, t) the quantum walk distribution 
and when all non-mirror paths (i.e. a ^ a') are non proper then 
p(s,t) — » PRw(Sjt) the classical random walk distribution. 

To probe the behaviour of p(s,t) we initially perform ex- 
act numerical simulations. Note that z(a,a!) is a sum of local 
characteristics of the links, while %(a,a!) is a sum of non-local 
ones. To compute quantities which involve pairs or triples 
of components, we make use of the braid word representa- 
tions of the links Eqs. (2,3) determined by the coin histo- 
ries {a, a!). The linking number lk(L Wl Lj) = (#bj — #bj)/2 
when the walker's component is to the right of j (w > j) 
and lk(L w ,Lj) = (#bj-#b])/2 otherwise. If lk{L w ,Lj) is an 

odd integer for any component Lj then L = (Bl,^ Bl) Markoy 
is not proper. The Milnor invariant C3(L rj L Wj L s ) where the 
components are ordered in position with r < w < s, can be 



4 





30 r 




25 ■ 




20 ■ 


CM 

b 






10 ■ 




5 ■ 




- 



• f 



• + 

X 



Ising 
RW 
QW 



10 15 20 25 
t 



FIG. 3: The time evolution of the variance of the spatial distribu- 
tion p(s,t) (Eq. 7) for the Ising anyonic walk and the corresponding 
classical and quantum walk evolutions with the same initial state. 



computed by writing down the braid word (Bl,^ Bl) Markov and 
deleting all instances of braid word generators not in the set 
{b r , b s - 1 , b\ , b\_ x } . The remaining braid word B (L r , , L s ) 
involves just two generators and their adjoints. We can think 
of this as a braid word on the three strands braid group. For the 
spin- 1/2 irrep of the Ising model, the representation of these 
generators is: 



1 

i 



V2 



1 



1 



Note that b 2 r = e - in l A G z and b]_ x 



e~ in / 4 a x where a xz are the 
Pauli operators. For proper links the product of this represen- 
tation of the braid generators in B(L ri L Wl L s ) will be ±1 with 
the sign carrying the value ( — l)^^' 1 *^). The Milnor invari- 
ant for the other two orderings of components is calculated 
analogously. We have calculated the distribution for walks up 
to t = 25. The variance c 2 (t) = (s 2 ) — (s) 2 where the expec- 
tation value is (0(s)) = £ 5 0(s)p(s,t), is plotted in Fig. 3 and 
it quickly approaches the linear random walk variance. Using 
the total variation distance between two distributions p(s,t) 
and f(s,t) denned A(p,f) = \p(s,t) -f(s,t)\, at t = 25 
we find A(/?,/?qw) = 0.34 while A(/?,/?rw) = 0.04. 

We now show that asymptotically the Ising anyonic walk 
behaves classically. The essential reason is the rapidly de- 
creasing density of proper links in the regime where the quan- 
tum walk distribution has dominant support. In order to up- 
per bound the variance we can assume that %(ci,d f ) of all the 
proper links are even and that there is no correlation between 
being proper and z(a : a f ) for non mirror paths. Calling the 
resulting distribution p(s,t)we have 

p(s,t) = p(s,t) RW + p prov (s,t)[p(s,t) QW - p(s,t) RW )}, (8) 

where Pp T op(s,t) is the density of proper links for non-mirror 
paths (a ^ a'). Since the walker's speed is constant, the maxi- 
mum possible variance is quadratic, achieved up to a constant 



less than one by the quantum walk [22], so this choice of dis- 
tribution can only make the estimate of the variance of the 
anyonic walker larger, i.e. d 2 (t) > o 2 (t). 

In Ref. [ ] it was shown that the distribution of the 
quantum walk with the same Hadamard coin flip and initial 
state as occurs here can be very well approximated asymptot- 
ically by the function p' n ^(aj)da = (i-a)^g — = where 

J J ^ QWV ' 7i(l-a 2 )x/l-2a 2 

a = (so — s)/t and is restricted to the interval [— 4j , 4j] . Out- 
side this interval the distribution falls off exponentially with 
t as does pqw(s,t). Restricting to this interval and using the 
fact that a proper link must have all linking numbers even, it is 
shown in the Appendix that the density of non-mirror proper 
links is Pprop(s,t) <C/t 2 for some constant C independent of 
s,t. The position moments with respect to p(s,t) are 

((s-So)) ~ f~ l ^p pmp (sj)tap f QW (aj)da 

< C(l-l/>/2)A, 

((s-S ) 2 ) w ^(l-CA 2 )+/-^ prop (^ 2 a% w (a,0^ 

< ?(l-CA 2 ) + (l-l/x/2)C, 

Thus d 2 (t) <t + 0(l). 

Similarly, to obtain a lower bound for the dispersion we can 
assume a distribution of the form of Eq. 8 but chose a prob- 
ability distribution f(s,t) with minimum variance to replace 
pqw(s,t). Picking f(s,t) = S s ^ (zero variance) and calling 
the resulting distribution p(s,t) then b(t) 2 > t(l-C/t 2 ). By 
the inequalities b 2 (t) < G 2 (t) < a 2 (t) we find 



lim 



a 2 (t) 



1, 



(9) 



i.e. asymptotically the Ising anyonic walk has linear disper- 
sion with coefficient 1 . 

In conclusion, we have studied the dynamical behavior of 
a mobile non-Abelian anyon which becomes entangled with 
its environment purely by statistical interactions. We find that 
for the case of Ising anyons the decoherence is strong enough 
to completely wash out the quantum mechanical interferences 
and reduce the dynamics to a classical stochastic process. It 
would be of interest to extend this analysis to anyons which 
are spin- 1/2 irreps for SU (2)^>2 models since it is known that 
for k > 2,k 7^ 4, the braiding evolutions densely span the fu- 
sion space while for k = 2 (Ising anyons) and k = 4 they do 
not [23]. It has been shown that for a t step quantum walk 
subject to decoherence in its position at a rate p mQas >C'/t 
for some constant C', the evolution approaches classical be- 
havior [ ]. We might then conjecture that since the braid- 
ing generators that entangle new fusion degrees of freedom 
can be translated into a measurement error p meas ~ l/k 2 (for 
k ^> 1), then for t > k 2 the walk would behave classically. 
In the limit k — >• °o the particles are fermions and we recover 
quantum walk behavior with quadratic dispersion. In general 
a mobile anyons experiences an environmental coupling that 
is highly non-local and non-Markovian and its behavior sheds 
light on the statistical dynamics of more complex systems. 
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Appendix: Upper bound for p vmv (s,t) 

We will show that p vmv (so : t), the probability of a randomly 
chosen link L to be proper (mirror paths excluded), is bounded 
from above by a function § under certain assumptions. We 
first consider the number of links with even linking num- 
bers and derive /? pr0 p 0*0,0 as the proportion of these links 
with respect to all links. The result holds for paths which 
end up on the initial site but numerical results suggest that 
Pprop(sjt) < /?prop(so,0 if s ^ so, provided s is not to close 
to the boundaries at so±t. In particular, it is true inside the 
domain *=* = [ --fe,-fe]. 

Suppose a path (3, a!) leading after t steps to the posi- 
tion so links with link components Lj, j G {/,..., r} (Fig.4). 
For the corresponding link L to be proper, all the linking 




FIG. 4: A generic path reaching after t steps (here t = 8) the position 
sq. The path braids with all the link components L/, ...,L r . We say 
that its width w is r — I + 1 (here w = 5). 

numbers lk(L w ,Li),...,lk(L w ,L r ) must be even. Let p Q (tJ) 
be the probability that lk(L Wl Lj) is even. We assume that 
for large enough (fixed) t, the probabilities p e (tj) can be 
treated as independent. This is justified inside the interval 
= [— -j=, -j=], where the number of paths that contribute 
to the anyonic walk density at each point s is exponential in t. 
The linking number of the walker with any particular compo- 
nent can change even/odd parity by a simple deformation and 
within the typical width (i.e. number of components touched 
by the walker) of the path there are an exponential number of 



such deformations hence the linking numbers of the walker 
with those components are well approximated as independent 
quantities. 

Then we can write the probability for the link L to be proper 

as 

/Vop0o,0 =Pe(t,l)'...'Pe(t,r) . (10) 

If we denote 

p = max p Q (t J) , (11) 
j 

where j runs over all link components Lj that braid with some 
t step long path leading to so, we can estimate from above 

PpropM<p W , (12) 

where w — r — / + 1 we will call the width of the link. 

The paths (3,3') relevant for a walk of t steps must satisfy 
a t = a' t . Let us consider the paths leading to so (we presume 
t is even), for which a t = a' t = 1. It is useful to depict 3 as a 
lattice path on a n — 1 x n lattice (n = |) with allowed steps t 
and (Fig. 5). The number of all such paths is 




FIG. 5: Any half -path (i.e. either "bra" or "ket" history of the 
walker) can be thought of as being realized on a two dimensional lat- 
tice, where only the steps t and — > are allowed. If we fix the number 
of steps t (presume t even), the position reached after t steps s = so 
and the last coin outcome a t = 1, we can draw the possible half-paths 
on a lattice n — 1 by n, where n = | . The lattice paths that touch the 
diagonal y = x + w (1 < w < n), but don't touch y = x + w + 1 , braid 
with at least w link components. "At least", because they certainly 
braid with the w components to the left of so, but they might be (and 
typically are) spread to the right of so as well. 



(^_^=#(all). (13) 

The number of lattice paths that touch the diagonal y = x + 
w, w G { 1 , . . . , n} is in our case of n — 1 by n lattice 

P^W). (14) 
\n — wj 

To derive this result we notice that #(1) = #(all) and for w G 
{2, ...,ft} we use the Andre reflection principle, described, for 
example, in [25]. 
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Consider the set fP w of lattice paths that touch the diagonal 
y = x J r w but do not touch the further diagonal y = x + w + 1 . 
The sets @ w are for distinct w G {1, ...,ft} obviously disjoint. 
The number of paths in !P W is 



|# w |=#(w)-#(w+l) 



2ft -1 
ft — w 



In -I 

n — w — 1 



(15) 



The number 2 in both the nominator and the denominator ac- 
counts for the fact that we count both types of paths - with 
a t = a' t = and a t = a! t = 1 - and use the symmetry between 
the two situations. We used the relation (17) and the diagram 
at Fig. 6. Also, we didn't forget to exclude the mirror paths. 
Before substituting to (18), we rewrite (15) as 



We use the convention that (™) = if k < 0. Furthermore, the 
union LX=i comprises all the paths on n — 1 by n lattice, 
which follows from 



(2?z- 1)! 



1 



1 



II*w|=I 



w= 1 



w= 1 



In -1 

ft — w 



2ft -1 
ft — w — 1 



2ft- 
ft - 



= #(aH) . 



(ft-w- l)!(ft + w- 1)! 
and (13) as 

(2ft- 1)! 



ft — w n-\-w 



w ( 2ft 



#(all) 



(16) 

A path (3,3') consists of two lattice paths, 3 and 3', where 
a e !P W and 3' G fP w / for some w, w' G {1, ...,ft}, i.e. (3,3') G 
x Let's realize that the width of the path (a, a 7 ) is at 
least max{w, w'} (see again Fig.4 and Fig. 5). The probability 
that the path (3, 3') is proper can thus be estimated from above 
(using (12)) as 



Ww)<p max{wX} 



(17) 



Finally, the set of all non-mirror paths with the last step "1" 
reaching after t steps the starting position so can be expressed 
in terms of fP w 's as shown in Fig. 6. 



-Pi 



mirror paths 
are excluded 



To 



FIG. 6: The set of paths that after t steps reach the posi- 
tion s = so, with the last coin outcome a t = 1, can be ex- 
pressed in terms of the sets !P W as [U w w > x ^V] W = 

[UU x <P W ) U (Up/ fPy x fP w ) U (u^ 1 x ft)] \M, 
where M = {(a,a')\a = a'}. 



Now we can write an upper bound for p vrov (so,t) as the 
number of links with even linking numbers divided by the 
number of all links: 



/VopC^OjO 



< 



-(g = iP w (i 



\tPw 



-2IJ- 



1 \- L w 



•a)) 



2(#(a//) 2 -#(all)) 



(18) 



ft!(ft-l)! 



1 /2ft 

2 V n 



n \n — w 
(19) 



(20) 



Now (18) reads 



Pprop(^O)f) < J737 



( 2 „") SUi p w (C-J jjgvl+zEgj W) 

2 ( 2 ;)-i 



1 /2n\ 1 /2n\ 
2UJ2UJ 



which we estimate from above using the relation between bi- 
nomial coefficients ( ) < ( 2 ^) to get: 



^prop(^0,0 < 



(2n\ V n n w ( -4-9 



1 /2rc\ 

2 Vw, 



1 



1 



1 f2/? 



1 



1 f 2n 
2 



1 

^ w— 1 



11 

2 2 



:Ep m 

w= 1 
n 

Lp 



w-i \ 



w w 3 



(21) 



vv= 1 



where we used the formula J 
(i.e. f > 4) 



iw(w — 1). For ft > 2 



1 (In 
2 



< 



2 



(2) 



Thus 



f>prop(*0,0 < <- n 2 



6 4" w 3 64~ w3 



(22) 



(23) 



VV= 1 



The infinite sum on the right hand side converges for all p G 
[0, 1) and one can find that 



w w 3 



p(l+4p + p 2 ) 



(1-P) 4 ' 

Hence we obtain the following upper bound (recall n - 
1 6-4-4p(l+4p + p 2 ) 



t z 

C 



(1-P) 4 



(24) 
■■t/2): 
(25) 
(26) 



where C doesn't depend on time (number of steps) t. 
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